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^ ■ Abstract 

■^ ■ This report summarizes the final results from the OPAL collaboration on searches for 

neutral Higgs bosons predicted by the Minimal Supersymmetric Standard Model (MSSM). 



^ . CP-conserving and, for the first time at LEP, CP-violating scenarios are studied. New 

' ! scenarios are also included, which aim to set the stage for Higgs searches at future colliders. 

(D \ The results are based on the data collected with the OPAL detector at e+e~ centre-of- 

i-^ ■ mass energies up to 209 GeV. The data are consistent with the prediction of the Standard 

Model with no Higgs boson produced. Model-independent limits are derived for the 

S^ . cross-sections of a number of event topologies motivated by predictions of the MSSM. 

Limits on Higgs boson masses and other MSSM parameters are obtained for a number of 

representative MSSM benchmark scenarios. For example, in the CP-conserving scenario 

?7ih— max where the MSSM parameters are adjusted to predict the largest range of values 

for rrih at each tan/3, and for a top quark mass of 174.3 GeV, the domain 0.7 < tan/3 < L9 

is excluded at the 95% confidence level and Higgs boson mass limits of ?7ih > 84.5 GeV and 

ruA > 85.0 GeV are obtained. For the CP-violating benchmark scenario CPX which, by 

construction, enhances the CP-violating effects in the Higgs sector, the domain tan/3 < 2.8 

is excluded but no universal limit can be set on the Higgs boson masses. 
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1 Introduction 

It is generally assumed that the Higgs mechanism p] is responsible for the breaking 
of electroweak symmetry and for the generation of elementary particle masses. In its 
simplest form, implemented in the Standard Model (SM), this mechanism implies the 
self-interactions of one doublet of complex scalar Higgs fields and predicts the existence 
of one physical scalar particle, the SM Higgs boson. The mass of this particle is not 
predicted by the model. Despite an intense effort over the last few decades, the SM Higgs 
boson has not been detected; a lower bound of 114.4 GeV has been set on its mass [2j. 

The SM, while successfully describing all electroweak phenomena investigated so far, 
provides no mechanism for stabilizing the electroweak energy scale (~ Mw) in the pres- 
ence of quantum corrections (the "scale hierarchy" problem). Supersymmetry (SUSY) is 
proposed as a possible solution to this problem since in SUSY models fermionic quantum 
corrections are compensated by bosonic corrections of similar size, and vice versa. 

The Minimal Supersymmetric Standard Model (MSSM) is the SUSY extension of the 
SM with minimal new particle content. It introduces two complex Higgs field doublets 
which is the minimal Higgs structure required to keep the theory free of anomalies and 
to give masses to all fermions. The MSSM predicts five Higgs bosons: three neutral and 
two charged particles. At least one of the neutral Higgs bosons is predicted to have its 
mass close to the electroweak energy scale; when radiative corrections are included, this 
mass should be less than 135 GeV |3j. This prediction provides a strong motivation for 
searches at present and future colliders. 

In the MSSM the Higgs potential is invariant under CP transformation at tree level. 
However, it is possible to break CP symmetry in the Higgs sector by radiative correc- 
tions, especially by contributions from third generation scalar-quarks |1]. Such a scenario 
is theoretically attractive since it provides a possible solution to the cosmic baryon asym- 
metry [5^ while the size of the CP-violating (CPV) effects occurring in the SM are far too 
small to account for it. 

The searches performed so far at LEP have been restricted to CP-conserving (CPC) 
MSSM scenarios where all CPV phases in the soft SUSY-breaking Lagrangian related 
to the Higgs sector are put to zero. Under this assumption the three neutral Higgs 
bosons are CP eigenstates: the h and H bosons (h is defined to be the lighter of the two) 
are CP-even and the A boson is CP-odd. Of these, only the CP-even states couple to 
the Z boson in Higgsstrahlung at tree level; thus, at LEP energies, these particles are 
mainly produced through the Higgsstrahlung processes e"'"e~— >^hZ and HZ, and the pair 
production processes e"'"e~— s>hA and HA (e"'"e~— s>hH is forbidden by angular momentum 
and CP conservation). Production of h and H via WW and ZZ fusion in the t-channel, 
where a Higgs particle is produced in association with a pair of neutrinos or electrons, 
is included and plays a role at the kinematic limit of Higgs boson production at LEP 
energies. 

In this paper CPV MSSM scenarios are also considered. In such scenarios the three 
neutral Higgs bosons, Hj {i = 1,2,3) are mixtures of the CP-even and CP-odd Higgs 
fields. Consequently, they all couple to the Z boson by their CP-even field components 
and to each other. Therefore, the Higgsstrahlung processes e"'"e~^H,jZ {i = 1,2,3) and 
pair production processes e"'"e~— i>HjHj {i 7^ j) may all occur, however with widely varying 
cross-sections. In large domains of the model parameters, the lightest neutral Higgs boson 
Hi may escape detection, since its coupling to the Z boson may be too weak. The other 
two Higgs boson masses may be out of reach on kinematic grounds or may also have small 



production cross-sections. As a result, the exclusion limits obtained previously for the 
CPC scenario may be invalidated by CPV effects. 

The decay properties of the Higgs bosons, while being quantitatively different in the 
two scenarios, maintain a certain similarity. Since Higgs bosons, in general, couple to mass, 
the largest branching ratios are those to bb and t^t~ pairs. If kinematically allowed, the 
cascade decays h— >AA (CPC scenario) [Sj or H2— >HiHi (CPV scenario) would occur and 
can even be the dominant decays. The searches described below take this possibility into 
account. 

Since the topological searches are essentially insensitive to the CP parity of the Higgs 
bosons, their results can be applied both to the CPC and CPV MSSM scenarios. For 
example, the search for the Higgsstrahlung signal topology e"'"e~^hZ in CPC scenarios 
can be interpreted as a search for the processes e"'"e~— s>HiZ and H2Z in CPV scenarios; 
similarly, the search for the pair production process e"'"e~— >^hA in CPC scenarios can be 
interpreted as a search for the process e"'"e~^HiH2 in CPV scenarios. 

The results presented in this paper are based on all data collected by the OPAL 
collaboration, including LEPl data taken at the Z resonance and LEP2 data collected 
between 130 and 209 GeV. Many of the searches are already described in earlier OPAL 
publications [7|IHll^ lTUl[mil21ll^ . In this paper we describe only the optimizations of the 
searches for the data collected at the highest energies, between ^/s = 192 and 209 GeV. 
The old and new searches are combined to produce the final results presented here. 

Recent searches for neutral Higgs bosons performed by other LEP collaborations, 
limited to the CPC MSSM scenarios, are listed in J3] while results from the Tevatron 
collaborations are reported in [T3] . 

This paper is organised as follows. Section|21summarises the main features of the OPAL 
detector and the Monte Carlo simulations. Section IHl describes the model predictions for 
Higgs boson production. The search topologies and the experimental searches for the 
corresponding final states are described in Section 01 with emphasis on those searches 
which have not been covered in earlier OPAL publications. In Section El the statistical 
combination of individual search channels and the derivation of exclusion limits in the 
MSSM parameter space are described. In Section IHl the search results are translated 
into model-independent limits on topological cross-sections. In Section (3 the MSSM 
benchmark scenarios, including those recently suggested in p!Uj, and the corresponding 
results are presented. A short summary is given in Section |H1 

2 OPAL detector and Monte Carlo samples 

The OPAL detector is described in jTTj. The tracking detectors and calorimeters have 
nearly complete solid angle coverage. The central tracking detector is placed in a uniform 
axial magnetic field of 0.435 T. The innermost part is occupied by a high-resolution 
silicon strip ( "microvertex" ) detector p^. It surrounds the beam pipe and covers the 
polar angle^ range | cos^| < 0.93, and is the basic tool for an efficient b-tagging p. This 
detector is followed by a high-precision vertex drift chamber, a large-volume jet chamber, 
and chambers to measure the z coordinates along the particle trajectories. A lead-glass 

^OPAL uses a right-handed coordinate system where the +z direction is along the electron beam and 
where +x points to the centre of the LEP ring. The polar angle 6 is defined with respect to the +z 
direction and the azimuthal angle with respect to the +x direction. The centre of the e'^e^ collision 
region defines the origin of the coordinate system. 



electromagnetic calorimeter with a presampler is located outside the magnet coil. The iron 
return- yoke of the magnet is instrumented with streamer tubes and thin-gap chambers for 
hadron calorimetry. Finally, the detector is completed by several layers of muon chambers. 

A variety of Monte Carlo samples are generated to estimate the detection efficiencies 
for the Higgs boson signal and to optimise the rejection of the background. The production 
cross-section and kinematic properties of the Higgs boson signal vary rapidly with energy 
near the kinematic limit. For an accurate modelling, the signal and background samples 
are generated at several centre-of-mass energies, from 192 GeV to 210 GeV. The main 
background processes are generated with 15 to 25 times the statistics of the data. 

Signal events for the Higgsstrahlung and pair production search channels are generated 
using the HZHA03 ^Hj program. The pair production samples are described in this paper. 
For Higgsstrahlung, the Monte Carlo samples used are described in [71. Typical sample 
sizes are 2000 to 5000 signal events per Higgs mass point. For the background processes the 
following event generators are used: KK2f pB^ for (Z/7)*^qq(7), fi~^fi~{'y) and r^T^{j), 
grc4f in] for four-fermion processes (4f), BHWIDE [22] for e+e-(7), and PHOJET |2S], 
HERWIG j23) s-iid Vermaseren J2S1 for hadronic and leptonic two-photon processes (77). 
JETSET [2^ is used as the principal model for fragmentation and hadronisation. The 
detector response is simulated in full detail fI7\ . 

3 Higgs boson production processes 

In the MSSM, the Higgsstrahlung and pair production processes have complementary 
cross-sections. Their relative rate is regulated by sum rules which are different in the 
CPC and CPV scenarios and depend on the precise choices of the MSSM parameters. 

In the CPC scenario, the cross-sections for the processes e"'"e~— >hZ, e"'"e~— s>HZ and 
e"'"e~^hA are given by 



e"'"e — i>hZ 
e+e-^HZ 
e+e"^hA 



ahz = sin^(/?-aXz(^h), (1) 

aHz = cos2(/5-a)4z('^H), (2) 

cThA = cos^(/3 - a) A 0-I2 (mh), (3) 



where a^ is the cross-section for the SM Higgsstrahlung process e"'"e~^HsMZ, tan/3 = 
^2/^1 is the ratio of the vacuum expectation values of the two Higgs field doublets coupling 
to "up" (^2) and "down" (vi) type fermions, and a is the mixing angle describing the 
combination of the two CP-even weak eigenstates to produce the two CP-even Higgs mass 
eigenstates. The symbol A denotes the kinematic phase-space factor 

x3/2 



A^/'(12M|/. + Azh) 

with 

_ 1 - (nii + rrijY Is 

^'' 1 — {rrii — rrijY / s 

Due to the complementarity of the Higgsstrahlung and pair production processes, ex- 
pressed in Equations (0) and Q, the searches have to include both of them in order to 
maintain a high sensitivity over the whole MSSM parameter space. 



Similar, but more complex, sum rules regulate the relative rates in the CPV scenario. 
The cross-sections for the processes e^e~— >-HjZ and e"'"e~— >HjHj are given by 



e+e-^H,Z : ^h.z = gl^zz 4z (^hJ, (4) 

e+e--^HiHj- : (Th.h, = ^h,h,z ^ ^Iz ("^hJ, (5) 



where the couphngs 



S'H.zz = COS 13 On + sin (3 Os, 
5'H,H,z = Osi (cos (3 O2J - sin /5 0^) - Osj (cos /? O2J - sin /3 On 



obey the sum rules 



7 . 5'h,zz 



Iv^ 



5'HfcZZ - 2 Z^ ^ijkgH^HjZ- 

The orthogonal matrix Ojj (i, j = 1, 2, 3) relating the weak eigenstates to the mass eigen- 
states has non-zero off-diagonal elements only in a CPV scenario. 

4 Signal topologies and experimental searches 

The topological searches are devised to detect Higgs boson production via the Higgs- 
strahlung and pair production processes. Many of these searches are described in earlier 
publications. Here we emphasize the most recent searches, optimized at the highest LEP 
energies and some modifications of the earlier searches to increase the sensitivity in parts 
of the MSSM parameter space which were not sufficiently covered before. The topolog- 
ical searches for the Higgsstrahlung process are listed in Tabled and those for the pair 
production process in Table El 

In the following, the notation e'^e~^7iZ and e'^e~—^7ii7i2 is used as a generic notation 
designating both the CPC processes e+e~— >hZ, e+e~^HZ and e"'"e~— s>hA, and the CPV 
processes e^e~— >HiZ, e^e~-^H2Z and e"'"e~— >HiH2. Since in e"'"e~ collisions the kinematic 
properties of the CPC and CPV signal processes are expected to be very similar, the 
topological searches described in the following can be applied in both scenarios. In both 
CPC and CPV models, the production process in Higgsstrahlung only contains couplings 
between the Z boson and a CP-even Higgs state, and the pair production always involves 
a CP-even and a CP-odd state. The production angle distributions are therefore expected 
to be the same. In CPV models, the Higgs mass eigenstate is a mixture of CP eigenstates 
and is able to change from a CP-even to a CP-odd state. Even in this case, the Higgs 
decay angle distributions will stay the same since both the CP-odd and the CP-even states 
are spin bosons. Small differences not measurable at LEP may arise from different spin 
correlations in the decay products of the Higgs boson. 

4.1 Searches for Higgsstrahlung processes 

(a) The most efficient search for the Higgsstrahlung processes is the one which is de- 
signed to search for the Standard Model Higgs boson. This search ffl is interpreted 



here as a generic search for the corresponding MSSM process e"'"e~— >7^Z. It takes 
advantage of the preferential decay of Higgs bosons into bb and t^t~ pairs and ad- 
dresses the following Z boson decays: Z -^ qq, i/z/, e+e~, /i+/i~, t^t^ . Moreover, this 
search is also sensitive to contributions to the signal from the W"*"W~ and ZZ fusion 
processes e"'"e~— >7iz/z/ and Tie^e", which may become important at the kinematic 
limit of the Higgsstrahlung process. 

(b) The Higgs cascade decay l-L2^''HiHi may play an important role in regions of the 
MSSM parameter space where it is kinematically accessible. In order to increase 
the sensitivity to cascade decays, the search described in [7j is adapted, in those 
parts which deal with the "four-jet" final state e"'"e~— >-(7i— >-bb)(Z— >qq) and the 
"missing energy" final state e"'"e~^(7i^bb)(Z^z/z/). These searches modified for 
e'^e~^(TC2^'HiTCi)Z are described below in Sections 14.1.11 and 14.1.21 

(c) The search for Higgs cascade decays is complemented by an earlier search for 
e'^e~^(H2^'HiT-Ci)Z [29 , which is specifically designed to be efficient in the domain 
mn, < 10 GeV. 

(d) For Higgs bosons produced in Higgsstrahlung e~^e~^HZ and decaying into 
particles other than b-quarks or r leptons, a flavour-independent search for 
e"'"e~— i>(7-^— i>hadrons) Z [HIHI!?T] is used. 

4.1.1 Modification of the search in the four-jet channel 

The search for the SM Higgs boson in the channel e"'"e~— >-HsmZ— >bbqq [7] is modified 
to be sensitive to the cascade decay T-i2^-'HiHi- The event selection is identical to the 
SM search, and thus the same candidate events are observed with the same expected 
background. The whole event is forced into four jets using the Durham jet finder |H2j . 
If Til is not too heavy, the two jets from Tii— >bb are often joined into one jet. The SM 
four-jet search is therefore also efficient for this decay and the expected signal rates from 
e"'"e"-^(7i2— >'bb)Z and e"'"e~^(7i2^'^^i'^^i^bbbb)Z can simply be added. The efficiencies 
for various combinations of {m-^^^m-}^^) are given in Table El The shape of the distribu- 
tion of the discriminating variable T> ^ , however, differs for the two decay modes. T> is 
a product of a mass independent and a mass dependent likelihood. Depending on m-j-i^ , 
the mass reconstruction is diluted by wrong jet pairings inside one jet, and thus the like- 
lihood distributions are broadened. The signal distribution of T> is therefore constructed 
at each point of the MSSM parameter space, taking into account the changing relative 
contributions from the two decays by first adding the relative contributions of 7^2 ^bb 
and 7^2— *^bbbb in the two likelihoods and then calculating the product. The systematic 
uncertainties are essentially the same as for the SM channel e"'"e~— >-HsmZ— >bbqq [7]. 

4.1.2 Modification of the search in the missing energy channel 

For the data taken at a/s = 199 to 209 GeV, the Artificial Neural Network (ANN) analysis 
for e"'"e~— >HsMZ^bbz/z/ [7^ is reoptimized for 100 < mT^j < 110 GeV and modified to be 
sensitive to ?i2^bb and 7^2— ^'^^I'^^i— ^bbbb decays simultaneously. In this mass range, 
the TC2^'HiHi decay is crucial especially in the CPV scenario. 

In the preselection, the event sample is split into two subsamples, according to the 
2-to-3 jet resolution parameter 1/32 of the Durham jet finder. Subsample A contains events 



with j/32 < 0.05 (7^2— >bb events and most of the T-C2—^'HiT-Ci—*hhhh with hght Hi) and 
subsample B events with 1/32 > 0.05 (most of the 7^2— ^'^^I'^^i— ^bbbb events with heavier 
Til). The separating value of 0.05 is chosen so that the efficiency for H2— *-bb events in 
subsample A is approximately the same as in the SM search. Separate neural networks, 
ANNa and ANNb are then trained for events belonging to the two subsamples. For 
ANNa the training is based on 7-^2— ^-HiTii— >-bbbb signal events with m-^-^ = 12 GeV and 
100 < mn^ < 110 GeV, while for ANNg m^ = 40 GeV and 100 < mn^ < 110 GeV is 
used. The analysis is slightly different from the SM missing energy analysis Tj. 

The preselection cuts 1 to 5 are identical to those described in . They are designed to 
remove accelerator-related backgrounds (such as beam-gas interactions and instrumental 
noise), dilepton final states, two-photon processes and radiative qq events, and to select 
events with a significant amount of missing energy. The following new preselection cuts 
are applied: 

6. The tracks and clusters in the event are grouped into jets using the Durham algo- 
rithm. Depending on j/32, the event is either grouped into two jets (1/32 < 0.05) in 
subsample A or into four jets (1/32 > 0.05) in subsample B. Each event in subsample 
B is required to have at least one track per jet. 

Additional requirements are imposed for subsample A (7/32 < 0.05): 

7. The acoplanarity angle (180° minus the angle between the two jets when projected 
into the xy plane) must be between 3° and 100° to reject Z/7— >-qq events with 
back-to-back jets. 

8. To reduce the background from W"'"W~ events, the event must not have an identified 
isolated lepton ^ 



The effects of the cuts on the data, the simulated signal and background samples are 
given in TablelU The 12 (11) variables used as inputs to ANNa (ANNb) are listed below. 
All variables are scaled to values between zero and one, and in some cases the logarithm of 
the variable is used, which gives less peaked distributions and is therefore better suited to 
an ANN analysis. The distributions of some of these input variables are shown in Figs. ^ 
and 121 The following 9 variables are common to both ANNa and ANNb. 

1. The scaled effective centre-of-mass energy jHS| Vs'/y/s, 

2. The scaled missing mass m^iss/y/s, 

3. The polar angle of the missing momentum vector | cos^missl, 

4. The b-tag likelihood output Bi jH] of the first (highest energy) jet, 

5. The b-tag likelihood output B2 of the second highest energy jet, 

6. The angle between the first jet and the missing momentum vector, 

ln(l -COsZ(ji,Pmiss)), 

7. The angle between the second jet and the missing momentum vector, 

COsZ(J2,J3miss), 

8. The x^ of the one-constraint kinematic fit, In(x^z)) where the missing mass mmiss 
is forced to mz. 



9. The scaled missing momentum Pmiss/v^- 
The following three variables are used in the construction of ANNa 

10. The polar angle of the thrust axis | cos6'thr|, 

11. The acoplanarity angle of the jets In (0acop), 

12. The logarithm of the energy difference between the two jets In \Ei — E2\. 
The following two variables are used in the construction of ANNb 

10. The b-tag likelihood output B^ of the third jet, 

11. The b-tag likelihood output B4 of the fourth jet. 

The Higgs mass is reconstructed using the di-jet invariant mass after the 1 constraint 
kinematic fit. The distributions of the ANNa and ANNb output variables are shown in 
Fig. El and the distributions of the reconstructed mass in Fig. |3 Candidate events are 
selected if ANNa > 0.5 for events in subsample A and ANNb > 0.5 in subsample B. 
The efficiencies for signal events for both 7-^2— '•bb and 7i2^'^^i'^^i^bbbb are determined 
for both selections. In each point of the MSSM parameter space, the expected signal 
distributions in the output variables are added from both signal sources according to 
their expected rates. The number of candidate events in subsample A(B) is 11(8) with 
10.0 (7.2) events expected from background (see Table E}. The signal efficiencies are 
shown in Table El for various values of {m-H^ , m-?^^ ) . The reconstructed masses and the 
ANN outputs are used to construct the discriminating variable V, which is used in the 
statistical combination with other search channels. 

The systematic uncertainties for this channel are evaluated by analogy to the SM 
missing energy search [7j. They amount to 2.2% (9.0%) in the signal (background) for 
subsample A and to 2.5% (19.3%) in subsample B. The strongest contributions arise 
from the uncertainty of the B hadron fragmentation and from Monte Carlo statistics. 
Systematic uncertainties of 9% for the background and 1% for the signal are added in 
quadrature to account for the uncertainty introduced when simulating the separation of 
the events into the two subsamples. These are estimated by shifting the value of 1/32 of each 
event by the difference in the mean values of the background and the data distributions 
of ?/32 and repeating the selection with the modified value of 1/32. 

4.2 Searches for pair production processes 

(a) The search for the four-b final state e"'"e~^(7ii^bb)(7i2^bb) provides the highest 
sensitivity. While in the CPC scenario the pair production process is dominant only 
for m-H-^ ~ m-Hi , this is not the case in the CPV scenario. The search in this channel 
is therefore optimized separately for small rrifi-^^ and large 7717.^^ . These are described 
below in Sections KTT\ and \rn\ 

(b) For the Higgs cascade decay e^e^^(Hi—*hh){H2~*'HiTi.i—^hhhh) with 6 b-quarks 
in the final state, the search in the four-b final state for similar masses described in 
Section B:.2.1l is used because it has a reasonably good efficiency. Even for large mass 
differences and thus small rriT-i^ this search is more efficient than the one described 
in Section I4.2.2[ due to the highly spherical shape of the six-b events. This search 
is described in Section 14.2.31 
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(c) The search for the final states e'^e~^(Hi^hh){H2^T~^T~) and 
e"'"e~—>-(?-^i— i>r+r~) (7-^2— >bb) follow the technique described in |7| for the 
corresponding Standard Model channels. The final likelihood selection and its 
optimization for the MSSM case is described in Section [4.2.41 

4.2.1 Search for e+e~— >?^i7i2— ^bbbb optimized for high m->^^ 

Events from the process 7ii7i2^bbbb with high m-Hi have four energetic b-jets and a total 
visible energy close to the centre-of-mass energy. The dominant backgrounds arise from 
the four-fermion processes e"'"e~— >-ZZ and e"'"e~— *W^W~ and from two-fermion processes 
e"'"e~^qq(7). The events are forced into four jets using the Durham jet finding algorithm 
and the following preselection is applied. 



1. The event must qualify as a multi-hadronic final state according to 

2. The effective centre-of-mass energy ys^ is required to be higher than 0.794 y^, 

3. The 3-to-4 jet resolution parameter 2/43 [32] is required to be larger than 0.003, 



4. The C-Parameter j34], which is a measure of the spherical shape of the event, is 
required to be larger than 0.45, 

5. The sum of the number of reconstructed tracks and electromagnetic clusters not 
associated to tracks |HS1 belonging to each jet has to be larger than six, 

6. To discriminate against poorly reconstructed events, a 4-constraint kinematic fit is 
applied, using energy and momentum conservation; this fit is required to converge 
and the x^ probability is required to be larger than 10~^. 

For the events passing the preselection, a likelihood function is constructed from seven 
input variables, allowing the events to be classified as signal or background. The likelihood 
variables are: 

1. The four b-tagging discriminants Bi [H| for each of the four jets, ordered by energy, 

2. The logarithm of the jet resolution parameter 2/43, 

3. The event thrust value T, 

4. The estimate of the TC1TC2 production angle, | cos^dijetl, which is defined as follows. 
For the jet pairing that yields the smallest difference between the two dijet-masses, 
I cos^dijetl is the absolute value of the cosine of the dijet polar angle. 

The signal reference histograms are obtained using Monte Carlo samples with m^^j > 
60 GeV and rrii^,^ > 60 GeV. For the background, reference histograms are formed from 
e"'"e~— >-qq(7) and e"'"e~— >-qqqq events. The distributions of these variables are shown in 
Figure El Events are selected if they satisfy £ > 0.95, which provides the best sensitivity 
measured in terms of s/\/b + 2 for m-^^ = m-^j = 90 GeV. 

The numbers of observed and expected background events after each preselection cut 
and the final likelihood cut are shown in Table IHl for data taken at ^/s = 192 to 209 GeV. 
The distribution of the likelihood output is shown in Figured (a). The efficiencies for 
various combinations of {mT-i^^m-H^) are shown in Tabled 
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The mass of the Higgs boson candidates is reconstructed using a constrained fit re- 
quiring energy and momentum conservation. Figures IHfa)-(c) show the distributions of 
the sum of the reconstructed Higgs boson masses, Mgum = f^^^t + to^-^ , for the jet combi- 
nation with the largest, second largest and smallest value for |m^^ ~^h2\- No significant 
excess over the expected background is observed. The discriminating variable "D is a 
two-dimensional array of reconstructed masses m:^^ + '^t?! ^^d m^-l^^ — mi^^. 

The systematic uncertainties on the signal efficiencies and background expectation 
for the 7ii7i2— >bbbb search are given in Table |H1 They are evaluated as for the SM 
searches in [Zj and include Monte Carlo statistics (uncorrelated between channels, energies 
and signal and background), detector modelling, such as tracking resolution in r0 and 
z, hit matching efficiency in the silicon microvertex detector for rcf) and z, B-hadron 
decay multiplicity and momentum spectrum, c-hadron momentum spectrum, comparison 
between different SM Monte Carlo generators, uncertainties of the four-fermion cross- 
section (all taken to be fully correlated between channels and energies) and uncertainties 
in the modelling of the likelihood variables (taken to be uncorrelated between channels 
and fully correlated between energies of the same channel). The systematic uncertainty 
amounts to 3.1% for the signal and 10.3% for the background expectation. 

4.2.2 Search for e+ e' ^Hx'H^^hhhh for low m-n. 

The region 12 GeV < rrifi^ < 30 GeV and rriT-i^ > 90 GeV is of particular interest in the 
CPV scenario. The following selection is optimized for that kinematic region and replaces 
the selection of Section 14.2.11 Events with large mji^ — m-j^^ look like asymmetrically 
boosted three-jet events. 

The preselection is identical to the one of Section f4.2.1l except for cuts (3) and (4). 
In (3) the cut on 1/43 is relaxed to 0.0003. In order to compensate for the increased 
qq7-background an additional requirement is introduced in cut (3): the sum of the two 
smallest angles between any jets, J2, has to satisfy the requirement 30° < J2 < 175° and 
the sum of the four smallest angles between jets, J4, has to satisfy 220° < J4 < 400°. 
Cut (4) is relaxed to C > 0.2. This increases the acceptance for asymmetric three-jet-like 
events. The number of selected events after each cut, along with the expected background, 
is shown in Table IHl 

After the preselection, a likelihood function is constructed from the seven variables 
described in Section 14.2.11 Signal reference histograms are formed from Monte Carlo 
samples with 12 < m-j-i^ < 30 GeV and 90 < m-)^^ < 110 GeV. The distributions of 
the input variables are shown in Fig. IHl The resulting likelihood distribution is shown 
in Fig. [7|^b). The cut C > 0.98 is applied, which is optimal for m-j-c^ = 30 GeV and 
m^2 = 100 GeV. 

The efficiencies for various combinations of {my^-^^niT-i^) are shown in Table ITUl The 
distribution of the reconstructed mass sum Msum is shown in Fig. IHl No significant excess 
over the background is observed. The discriminating variable "D is a two-dimensional array 
of reconstructed masses m^^^+^Wi ^^^ ^^t^j^^^Wi- '^^^ systematic uncertainties are listed 
in Table |H1 and are derived in the same way as for the search described in Section 14.2.11 
They amount to 4.7% for the signal and 10.5% for the background expectation. 

4.2.3 Search for e+e-^TiiTis^Wi'HiTii^bbbbbb 

The search channel for 7Yi7Y2— >^bbbb optimized for high mji^ is also used to search for 
events of the type e+e~— >?-^i7-^2^'^i'^i'^i^bbbbbb. Despite the large mass difference 
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m-^2 ~ '"^Wi ^^d generally relatively low m-^-^, the selection for high 7717^^ f Section I4.2.ip 
is more efficient than the selection for low 7717^^ f Section I4.2.2|) due to the spherical shape 
of the bbbbbb events. 

The expected signal distribution is added to the one from the 4b-channel in the same 
way as described in Section 14.1. II The efficiency of this search is shown in Table ^2 for 
various [m-H-^jm-H^)- The systematic error of the signal of this channel is taken to be the 
same as for the search optimized for high m-^^. 

4.2.4 Search for the process e+e~— >?ii7i2 — >hh t~^t~ and r+r^bb 

The search for the 7ii7i2-^bbr+r^ final state, where either Hi or 0.2 decays into a tau 
pair, uses the same techniques as the SM tau search |7;. The final likelihood selection 
including the b-tagging is optimized for the MSSM process. The likelihood C'HiH2 differs 
from the SM likelihood in the following way. 

1. The signal reference histograms are constructed from simulated 7i\H2 events, using 
a large range of m-j-i^ and m-}^^ values, 

2. The x^ probability of the 3-constraint ffi, constraining the invariant mass of the 
T^T^ pair to mz, is dropped, 

3. The variable |cos6'dijet| (cf. Section l4.2.1|) is introduced. Here the dijet pairing is 
defined by the pair of b-tagged jets and the pair of tau jets. The mean value of the 
cosines of the two systems is taken. 

The distributions of the likelihood input variables are shown in Fig.^1 The distribution of 
the likelihood C is shown in Fig. ^2 ^ind the number of selected events after each cut along 
with the expected background is shown in Table IHl The discriminating variable D is a two- 
dimensional array of reconstructed masses m^-^^ +771^5^^ and nf-^^ — m^^^. The cut C > 0.64 
is applied. The efficiencies of the selection are given in Table IT^ Systematic uncertainties 
are listed in Table |H1 They are derived in the same way as in Section 14.2.11 Additionally 
the systematic error on the tau identification is evaluated using event mixing [2j. They 
amount to 2.6% for the signal and 15.8% for the background expectation. 

4.3 Comparison of Data and Expected Background Monte Carlo 

In Table El the results of the new searches, which are described in detail in this section, 
are summarized in terms of total background and data events. There is good agreement 
between the selected data events and the expected SM background. The numbers in this 
table are obtained by choosing for illustration my^-^ = 39 GeV and m-j^^ = 105 GeV in 
those searches where the selections depend explicitly on the hypothetical Higgs boson 
masses. 

No significant excess of candidate events over the expected background is found in any 
of the old and new search channels. 
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5 Combination of search channels and hypothesis 
testing 

5.1 Combined confidence levels 

The sensitivity of the searches for hypothetical Higgs bosons is increased by combining the 
results of the various topological searches. This is done following the statistical method 
described in ^7j. 

In order to compute confidence levels, a test statistic Q is defined that can be used to 
quantify the compatibility of the data with two hypotheses: the background hypothesis 
and the signal+background hypothesis. The confidence levels are computed from a com- 
parison of the observed test statistic and its probability distributions for a large number 
of simulated experiments for these two hypotheses. In this paper, the ratio Q = Cg+b/'^b 
of the likelihoods for the two hypotheses is chosen as the test statistic. The results of 
all search channels are expressed in fine bins of discriminating variables V, as defined 
in the descriptions of the individual searches. For each bin i oi T> three numbers are 
calculated: Sj, the number of expected signal events for a given set of model parameters 
(Higgs masses etc.), ftj, the number of expected background events, and n^, the number 
of observed events. Each bin is considered to be a statistically independent counting 
experiment obeying Poisson statistics. The test statistic can then be computed [7j as 

-2 In Q = 2 ^ Si - 2 ^ n, ln(l + s./h). 

i i 

The confidence level for the background hypothesis, CLfe, is defined as the probability 
to obtain values of Q no larger than the observed value Qohs-, given a large number of 
hypothetical experiments with background processes only, 

CLfo = P{Q < Qobs I background). 

Similarly, the confidence level for the signal+background hypothesis, CLs+bi is defined 
as the probability to obtain values of Q smaller than observed, given a large number of 
hypothetical experiments with signal and background processes, 

CLs+b = P{Q < <5obs|signal + background). 

In principle, CLs+b can be used to exclude the signal+background hypothesis, given a 
model for Higgs boson production. However, this procedure may lead to the undesired 
possibility that a downward fluctuation of the background would allow hypotheses to be 
excluded for which the experiment has no sensitivity due to the small expected signal 
rate. This problem is avoided by introducing the ratio 

CLs = CLg+b/CLb. 

Since CL;, is a positive number less than one, CL^ will always be greater than CLs+b and 
the limit obtained in this way will therefore be conservative. We adopt this quantity for 
setting exclusion limits and consider a hypothesis to be excluded at the 95 % confidence 
level if the corresponding value of CL^ is less than 0.05. 

The expected confidence levels are obtained by replacing the observed data configu- 
ration by a large number of simulated event configurations for the two hypotheses back- 
ground only or signal+background. These can be used to estimate the expected sensitivity 
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of a search and to compare the observed exclusion with the one expected with no signal 
present. 

The effect of systematic uncertainties of the individual channels is calculated using 
a Monte Carlo technique. The signal and background estimations are varied within the 
bounds of the systematic uncertainties, assuming Gaussian distributions of the uncer- 
tainties. Correlations are taken into account. These variations are convoluted with the 
Poisson statistical variations of the assumed signal and background rates in the confidence 
level calculation. The effect of systematic uncertainties on the exclusion limits turns out 
to be generally small. 

In case of overlapping channels, i.e. channels sharing a fraction of events, the approach 
described above is modified. Such a situation occurs for example for the Higgsstrahlung 
searches with and without b-tagging '7','HT^ or in the case of the Higgsstrahlung four-jet 
channel and the pair production four-b channel. We calculate the expected CL^ for each 
of the overlapping channels in turn, and retain only the channel that yields the smaller 
expected CL^. This procedure is repeated for each signal hypothesis. For different Higgs 
boson masses therefore different search channels give the exclusion. 

The same procedure is applied if two signal processes, for example TiiL and 7^22, can 
contribute to the same event topology, but at different mass values. In the case of the 
four-jet channel the selection procedure and discriminant variable T) depend on the Higgs 
mass hypothesis (test mass). Two different test masses have not only different signal dis- 
tributions Si but also different background and data distributions hi and n^. The selected 
events in searches for Ti^L and for 7^2 Z cannot be combined since the inconsistent back- 
ground and data distributions for the two hypotheses in general contain an overlapping 
sample of data events. Therefore only the hypothesis that yields the lower expected CL^ 
is retained. 

5.2 Additional experimental constraints 

If a given model for Higgs production is not excluded by using the search channels de- 
scribed above, the following additional constraints are considered: 

(a) The constraint from the measured Z boson decay width Fz: a model is regarded as 
excluded if the condition 

AF 
y aw,z(91.4 GeV) + V an^n, (91.4 GeV) > (Tz(91.4 GeV) — ^ 

^ J- z 

is satisfied using results from [SEl- The nominal LEPl centre-of-mass energy of 
sfs = 91.4 GeV is used. Fz is the total Z width, az is the total Z cross-section 
and AFz = 6.5 MeV is the maximum additional width that is compatible with the 
measured width, given the SM hypothesis (obtained from ZFITTER [^) at the 
95% CL. 

(b) The constraint from the decay mode independent search for e^e~— *^HZ jHEj: a model 
is regarded as excluded if 

an,z > Hi^hJ cr|z ("^wJ with mn, = rrin 

is fulfilled, where klm-H^) is the smallest scale factor for the SM Higgs production 
cross-section that is excluded at the 95% CL by this search. This criterion is used 
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for Tii = Hi and Hi = H2 and at y/s = 91.4, 183 and 206 GeV. The use of Z width 
constraints and decay mode independent analyses is especially helpful for the range 
rriy^-^ < 6 GeV. In the CPV scenarios all points excluded by the Z-width constraint 
turn out to be also excluded by the decay mode independent search. 

(c) The constraint from a search for Yukawa production of a light Higgs boson \^^: a 
model is regarded as excluded if the predicted value of the Yukawa enhancement 
factor ^ for Hi, multiplied with the branching fraction BR(7^i— i>r+r~), is larger 
than the smallest value excluded in PS^. In the case of the CPV scan, where Hi 
is composed of CP-odd and CP-even parts, the weaker of the two limits calculated 
for the Yukawa production of a CP-even or a CP-odd Higgs boson is used in the 
comparison. For CP-even Higgs bosons, ^ = — sin a/ cos j3, while for CP-odd Higgs 
bosons C, = tan/3 holds. This constraint is helpful in excluding models with large 
tan/?, 2'mj. < m-}^^ < 2mb and vanishing e"'"e~^7iiZ cross-section. 

(d) As an overlay, the constraint from measurements of inclusive decays of a b quark 
into an s quark and a photon BR(b^s7) is shown. A model is shown as excluded 
by this constraint if the corresponding branching ratio, calculated using ^Hj, falls 
outside the bounds 2.33 x 10"^ < BR(b^s7) < 4.15 x 10"^ (95% CL) [H]. This 
limit is used to constrain the CPC scenarios. 

6 Model-independent limits on topological cross- 
sections 

For the model-independent interpretation of the OPAL Higgs searches the scaling factor 



S95 



O'rei 



is computed, where CTmax is the largest production cross-section allowed at 95 % CL and 
(Jref is a reference cross-section. For Higgsstrahlung the SM cross-section asm is used as 
(Jref; for pair production the cross-section of equation (jHl) with cos^(/5 — a) = 1 is used. 
Initial-state radiation is included according to jl2]. 

Cross-section limits on the SM-like production and decay can be found in [7j and for 
flavour independent 7i— >qq decays in [31 . 

Fig. 1121 (a) shows S95 for the production process e"'"e^^7i2Z^7Yi7iiZ^bbbbZ. 
BR(H2^HiHi) = 1 and BR(?ii^bb) = 1 is assumed. The observed borders and dis- 
continuities stem from a number of different searches contributing and being sensitive in 
different mass ranges. For m-j^^ < 80 GeV, specific searches for this final state at 183 GeV 
provide a strong exclusion. For 80 < m-^^g < 100 GeV, only the Z— i>qq final state is used, 
giving a weaker exclusion. For 100 < m-j-i,^ < 110 GeV, the Z^uu final state is also 
employed. The limits are calculated for ttit^^ > 10.5 GeV only where the decay Hi—>-hh 
becomes kinematically possible. 

Fig. 1121 (b) shows S95 for the process e^e^-^HiH2~^hhhh. BR(7ii— >bb) = 
BR(H2— >bb) = 1 is assumed. The kinematic limit for ^/s = 206 GeV is indicated as 
a dashed line. Most searches apply only for m-}^-^ > 30 GeV. Below m-H^ = 30 GeV, only 
searches for pair production at ^/s = 183 GeV or lower contribute. Additionally, the area 
of m-Hj > 12 GeV and 90 < m-H^ < 110 GeV is studied in the data at ^i = 199 to 
209 GeV. 
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In Fig.^l(a) S95 for the process e+e —^7ii7i2-^hhT^ r is shown. The branching ratios 
are set to BR(7^i^bb) = BR(7^2^bb) = 0.5 and BR(7^i^r+r-) = BR{n2^T+T-) = 
0.5. The kinematic hmit for ^/s = 206 GeV is indicated as a dashed hne. The domain 
below TTiTYi = 30 GeV is covered only by data collected at ^/s = 183 GeV or lower. 

Fig. El (b) shows the exclusion region for the process 
e+e--^nin2^ninini^hhhhhh. BR(7^2^HiHi) = 1 and BR(7^i-^bb) = 1 are 
assumed. At 777-7^2 < §0 GeV the exclusion is stronger than for higher m-H^ due to 
dedicated searches at ^/s up to 189 GeV. Above m-^j = 80 GeV, only searches using data 
recorded with i/i of 199 to 209 GeV are available. 

7 Interpretation of the search results in the MSSM 

The presence of neutral Higgs bosons is tested in a constrained MSSM with seven param- 
eters. Two of these parameters are sufficient to describe the Higgs sector at tree level. A 
convenient choice is tan /3 (the ratio of the vacuum expectation values of the Higgs fields) 
and one Higgs mass; mx is chosen in the case of the CPC scenario and mH± in the CPV 
scenario. Additional parameters appear at the level of radiative corrections; these are: 
^susY, M2, /i, A, and rrig. All soft SUSY-breaking parameters in the sfermion sector are 
set to msusY at the electroweak scale. M2 is the SU(2) gaugino mass parameter at the 
electroweak scale and Mi, the U(l) gaugino mass parameter, is derived from M2 using the 
GUT relation Mi = M2(5sin^6'vi//3 cos^6'vi/), where 9^ is the weak mixing angle^. The 
supersymmetric Higgs mass parameter is denoted /x. The parameter A = At = A^^ is the 
common trilinear Higgs-squark coupling for up-type and down-type squarks. The stop 
and sbottom mixing parameters are defined as Xt = At — /xcot/3 and Xb = Ab — /itan/5. 
The parameter mg is the gluino mass. For the CPV scenario the complex phases related 
to At^b and mg are additional parameters. The phase related to At^b enters at one-loop 
level while the one related to mg enters as a second-order correction to stop and sbottom 
loops. Large radiative corrections to the predicted mass m^-j-i^ arise from scalar top loops, 
while the contributions from scalar bottom loops are smaller. 

The precise mass of the top quark has a strong impact on rriT^^; it is taken to be 
"^top = 174.3 GeV, the current average of the Tevatron measurements [43]- To account 
for the current experimental uncertainty, all MSSM interpretations are also done for 
"^top = 169 GeV and mtop = 179 GeV. 

Rather than varying all of the above MSSM parameters independently, we consider 
only a certain number of "benchmark sets" where the tree level parameters tan f3 and m,A 
(CPC scenario) or mH± (CPV scenario) are scanned while all other parameters are fixed. 
Results are presented for eight benchmark sets PHITE] in the CPC scenario and nine in the 
CPV scenario J2H] • Each scan point within a given benchmark set defines an independent 
realization of the MSSM (a model), which is tested by comparing its predicted observables 
(masses, cross-sections and decay branching ratios) with the experimental data. The 
parameters of the scans are summarized in Table ^^ 

For a given scan point the observables in the Higgs sector are calculated using two 
theoretical approaches. The FEYNHIGGS program [131 H^ is based on a two-loop dia- 
grammatic approach [TTI BS] and uses the OS renormalization scheme, while SUBHPOLE 
and its CPV variant CPH J2Hj are based on a one-loop renormalization group improved 



^Ms, M2 and Mi are the mass parameters associated with the SU(3), SU(2) and U(l) subgroups of 
the Standard Model. M3 enters only via loop corrections sensitive to the gluino mass. 
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calculation |iniEDlEIlE2] and uses the MS scheme. Both calculations give consistent re- 
sults although small differences naturally exist. Numerical values for parameters in this 
paper are given in the MS scheme. 

In the CPC case, the FEYNHIGGS calculation is retained for the presentation of the 
results since it yields slightly more conservative results (the theoretically allowed param- 
eter space is wider) than SUBHPOLE. Also, FEYNHIGGS is preferred on theoretical 
grounds since its radiative corrections are more detailed than those of SUBHPOLE. 

In the CPV case, neither of the two existing calculations is preferred a priori on 
theoretical grounds. While FEYNHIGGS contains more advanced one-loop corrections, 
CPH is more precise at the two-loop level. We therefore opted for a solution where, 
in each scan point, the calculation yielding the more conservative result (less significant 
exclusion) is retained. For illustration, the results from FEYNHIGGS and CPH are also 
shown separately for the main CPV scenario CPX (see Section 17^ . 

The limits obtained for the different benchmark sets are summarized in Table El 

7.1 CPC benchmark scenarios 

Of the eight CPC benchmark sets examined in this paper, sets 1, 3 and 6 have been used 
in the past. Scenarios 4 and 5 are motivated by experimental constraints on the branch- 
ing ratio of the inclusive decay b^s7 and recent measurements of the muon anomalous 
magnetic moment {g — 2)^. Benchmark sets 2, 7 and 8 are motivated by the fact that 
the planned Higgs searches at the LHC may have low sensitivity to detect Higgs bosons 
in these situations. The choice of parameters is summarized in Table E| 

In most cases, tan/5 is scanned between 0.4 and 40. For values below 0.4 the theoretical 
predictions become unreliable; for tan/3 larger than 40 the decay width of the Higgs 
bosons may become comparable to or larger than the experimental mass resolution, and 
the modelling of the signal efficiencies may loose precision. The value of m^ is scanned 
between and 1000 GeV. For values of m\ < 2 GeV, the branching ratios of the A 
become dominated by resonances and their calculation is unstable. However this area can 
be probed using direct searches for the heavier h boson, decay independent searches and 
Fz constraints. 

In general there is good agreement between the data and the background estimation, 
therefore limits on the MSSM parameters can be derived. In the CPC scenarios, the 
largest observed excess of the data over the background appears at m^ = 95 GeV in the 
mh— max benchmark set. At this point the excess is (1 — CLb) = 2.6 x 10^'^, corresponding 
to a significance of 2.8 a. It should be noted, however, that there is a large statistical 
probability of such an excess to appear somewhere in the parameter space under study. 

1. In the no mixing benchmark set the stop mixing parameter Xt is put at zero. The 
other parameters are fixed at the following values: mgusY = 1 TeV, M2 = 200 GeV, 
fi = —200 GeV. The gluino mass mg has little effect on the phenomenology of this 
scenario; its value is set to 800 GeV. 

The corresponding exclusion plots are shown in Fig. El The unexcluded region 
with 64 < TTih < 88 GeV and rriA < 43 GeV is due to the dominance of the cascade 
decay h^AA for which the search sensitivity is lower than for the h^bb and t^t^ 
channels. One should note, however, that in this domain the charged Higgs boson 
mass mii± is predicted to smaller than 81 GeV. This area is probed by charged Higgs 
boson searches fSS^, which will be further extended in the future (see Fig. [ill (d)). 
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The region with rr^h > 83 GeV and mx > 82 GeV is still unexcluded. In this domain, 
either the cross-section for Higgsstrahlung e"'"e~— >-hZ is small (sin^ (/? — «) is close to 
0, see Eq. [T]in Section EJ or the pair production process e+e"— >^hA is kinematically 
forbidden. 

Values of tan/? are excluded from 0.8 to 6.2. However, the tan/3 limit is strongly 
dependent on the top quark mass which was taken to be r/itop = 174.3 GeV. For 
^top = 179 GeV, the tan/3 exclusion is reduced to 0.8 < tan/3 < 4.7. If one 
disregards the unexcluded domain at low tan/?, the following lower bounds are 
obtained at the 95% confidence level: mh > 83 GeV and mp^ > 82 GeV. 

The constraint from the measured value of BR(b^s7) (see Section I372|) is indicated 
in Fig. El (b). 

2. The no mixing (2 TeV) benchmark set differs from the no mixing scenario in the 
flipped sign of /i (which is preferred by the current results on {g — 2)^) and by a 
larger SUSY mass scale Msusy = 2 TeV. The value of tan /3 is scanned only from 
0.7 to 40 due to numerical instabilities in the diagonalisation of the mass matrix for 
very low tan/?. Therefore the largest part of the unexcluded region of the no mixing 
case at low tan/3 is not probed in this scenario. 

The corresponding exclusion region is shown in Fig. El For f^A > 2 GeV, i.e. 
above the region of resonant Higgs boson decays, absolute limits can be set for the 
Higgs boson masses and on tan/3, which are rrih > 83.3 GeV, rriA > 84.3 GeV and 
tan/? > 4.2. If the unexcluded area at ttia < 2 GeV is also regarded, the exclusion 
in tan/? is 0.9 < tan/? < 4.2. The reduced tan/? exclusion with respect to the no 
mixing case reflects the increased value of Msusy- This limit is further weakened 
to tan/? > 3.2 for mtop = 179 GeV. 

The measurements of BR(b— >-s7) exclude the no mixing (2 TeV) scenario for mp^ < 
450 GeV, as can be seen in Fig. HSl (b). 

3. The mh— max benchmark set is designed to yield the largest range of mh for a given 
(mA,tan/3). This scenario is therefore the most conservative in terms of exclusion 
in tan/?. The other parameters are fixed as in the no mixing scenario, with the 
exception of the stop mixing parameter Xt = v6 TeV. 

The exclusion plots for this benchmark set are shown in Fig. El The following 
absolute limits are obtained at the 95% confidence level: mh > 84.5 GeV and 
niA > 85.0 GeV. Furthermore, values of tan/? between 0.7 and 1.9 are excluded. 
For mtop = 179 GeV this exclusion shrinks to the domain 1.0 < tan/3 < 1.3. Since 
the TTih— max benchmark set yields the most conservative exclusion in tan/3, also 
'Titop = 183 GeV as the anticipated 1 a upper bound of an increased world average of 
"^top = 179 GeV was tested. This is illustrated in FigEl(b) and in FigEl(c), where 
the exclusion in the (tan/3, mA) plane respectively the theoretical upper bounds on 
mh for mtop = 179 GeV and mtop = 183 GeV are also shown. Should the world 
average of the top quark mass move beyond 179.5 GeV, the exclusion in tan/3 would 
vanish completely. 

The supplementary constraint from the measured value of the BR(b— >S7) is shown 
in Fig. El (b). 
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4. The ?7ih— max"*" benchmark set differs from the ?7ih— max case only by the flipped 
sign of fi. This choice is favored by the presently available results on {g — 2)^ jMj . 
Since the Higgs boson properties depend only weakly on the sign of /i, the accessible 
Higgs mass range as well as the excluded domains are very similar to those of the 
mh— max scenario; they are shown in Fig. IT71 

The limits on the Higgs masses are rrih > 84.5 GeV and rriA > 84.0 GeV. The 
excluded range in tan/5 is 0.7 < tan/5 < 1.9, which decreases to 0.96 < tan/3 < 1.4 
for mtop = 179 GeV. 

The mh— max"*" scenario is excluded for rriA < 600 GeV by BR(b— >S7) measurements 
for all values of tan/? considered (between 0.4 and 40). This means that only the 
decoupling limit with mh at its maximum value for a given tan /3 is still allowed by 
BR(b^s7). 

5. The constrained m/^-max benchmark set differs from the mh— max"^ set by the flipped 
sign of Xt, which yields better agreement with BR(b -^ S7) constraints. One ob- 
serves that the maximum value of the Higgs boson mass at a given tan (3 is lowered 
by about 5 GeV. 

The excluded areas for this scenario (see Fig. fTH|l show similar features as the 
mh— max and ?7ih— max"*" scenarios. The limits on the Higgs masses are rrih > 
84.0 GeV and mp^ > 85.0 GeV. The excluded range in tan/? is 0.6 < tan/? < 2.2, 
which shrinks to 0.8 < tan/3 < 1.8 for mtop = 179 GeV. This is illustrated in 
Fig.UHl) (b) and (c). 

The supplementary constraint from the measured value of the b— >S7 branching ratio 
is shown in Fig. ^1 (b) as the band delimited by the two dash-dotted lines. 

6. The large fi benchmark set is designed to illustrate choices of parameters for which 
the detection of the Higgs bosons is believed to be a priori difficult at LEP. The 
parameters are set to the following values: m,susY=400 GeV, /i = 1 TeV, M2 = 
400 GeV, mg = 200 GeV, Xt = -300 GeV. It is scanned from tan/5 = 0.7-40 and 
mA = - 400. 

For this set of parameters, the h boson is always kinematically accessible (mh < 
108 GeV) but its decay to bb, on which most of the searches are based, is suppressed. 
For many of the scan points the decay h— i>r+r~ is also suppressed. The dominant 
decay modes are thus h^cc,gg or W+W~, and the detection of Higgs bosons has 
to rely more heavily on flavour-independent searches. 

In some of the scan points the Higgsstrahlung process e+e^— i>hZ is suppressed all- 
together (sin^ (/?-«) small). However, the heavy neutral scalar is relatively light in 
such cases (mn < 109 GeV) and the cross-section for the process e"'"e~— >HZ, being 
proportional to cos^(/5 — a), is large. 

The exclusions for this benchmark scenario are given in Fig. ^1 They show that 
the parameter space is essentially excluded even in this difficult scenario, with the 
exception of a few isolated "islands" . Those may slightly increase for higher values of 
the top quark mass. The origin of the islands can best be explained using Fig.ll9lfb). 
The large diagonal island at tjia > 100 GeV is due to the fact that BR(h^-bb) goes 
to there. The two thin vertical islands around ttia > 100 GeV are due to an overlap 
between e"'"e~— >^hZ and e^e~— i>HZ production. Both are kinematically accessible, 
but either one or the other can be used in the interpretation. 
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The supplementary constraint from the measured value of the b— s>s7 branching ratio 
is shown in Fig. ^J (b). 

7. The gluophobic benchmark set is constructed such that the Higgs coupling to gluons 
is suppressed due to a cancellation between the top and the stop loops at the hgg 
vertex. Since at the LHC the searches will rely heavily on the production of the 
Higgs boson by gluon-gluon fusion, such a scenario may be difficult to investigate 
there. The parameters chosen are : msusY= 350 GeV, M2 = 300 GeV, /i = 300 GeV, 
Xt = -750 GeV, 0.4 < tan/5 < 40, GeV < ttia < 1 TeV and mg=500 GeV. 

The exclusion for this benchmark set is shown in Fig. 1201 It is excluded to a large 
extend. The limits on the Higgs masses are rrih > 82 GeV and tua > 87.5 GeV. The 
excluded range in tan/3 is tan/? < 6.0. The excluded range is reduced to tan/3 < 3.5 
for rritop = 179 GeV. 

The supplementary constraint from the measured value of the b^s7 branching ratio 
is shown in Fig. 1201 (b). 

8. In the small Oefr benchmark set the Higgs boson decay channels h — i> bb and h -^ 
T^T~ are suppressed with respect to their Standard Model coupling by the additional 
factor — sinaefr/cos/3, stemming from corrections from b — g loops. This scenario 
may also be difficult to investigate by the LHC experiments. Similarly to the large- 
/i scenario, such suppressions occur for large tan/3 and not too large rriA- The 
parameters chosen are: msusY = 800 GeV, M2 = 500 GeV, ^ = 2 TeV, Xt = 
-1100 GeV, 0.4 < tan/3 < 40, and mg =500 GeV. 

The exclusion for this benchmark set is shown in Fig. j^ The limits on the Higgs 
masses are mh > 79.0 GeV and rriA > 90.0 GeV. The excluded range in tan/? is 
0.4 < tan/3 < 3.6, which is reduced to 0.5 < tan/? < 2.9 for rritop = 179 GeV. It 
appears that effects of suppression of the decays h^bb and h— >r+r^ do not play a 
role in the region kinematically accessible at LEP. 

The constraint from the measured value of BR(b— >S7) is shown in Fig. j^ (b). 

7.2 CPV benchmark scenarios 

In the MSSM, the Higgs potential is invariant under CP transformations at tree level. 
However it is possible to explicitly or spontaneously break CP symmetry by radiative 
corrections jHSI- In particular, the phases of At^h and mg introduce CP violation into the 
Higgs potential via loop effects, leading to sizeable off-diagonal contributions to the Higgs 
boson mass matrix (see Sectional). 

As a consequence the Higgs mass eigenstates Hi,H2 and H3 are not the CP eigen- 
states. This influences predominantly the couplings in the Higgs sector. Fig. (221 shows 
the coupling of a mixed mass eigenstate Hi consisting of admixtures from the CP eigen- 
states h, H and A. Since only the CP-even field component couples to the Z boson, the 
individual coupling of the mass eigenstates are reduced in the CPV case with respect to 
a CPC case. 

The size of the CPV off-diagonal elements, J^h, of the Higgs boson mass matrix and 
hence the size of CPV effect scales qualitatively |2^ as 



2 mf ImifiAt} 
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Large CPV effects, and thus scenarios dissimilar from the CPC case, are therefore obtained 
if the SUSY breaking scale frisusY is small and the imaginary contribution to ^At large. 
Also large values of mtop increase the CPV effects. 

When choosing the parameters, experimental constraints [SEIEZI from electric dipole 
moment (EDM) measurements of the neutron and the electron have to be fulfilled. How- 
ever, cancellations among different contributions to the EDM may naturally emerge [KH] : 
hence those measurements provide no universal exclusion in the MSSM parameter space, 
while direct searches at LEP provide a good testing ground for a CPV MSSM. 

The basic CPV MSSM benchmark set is CPX. Its parameters are chosen such as 
to approximately fulfill the EDM constraints and to provide features that are the most 
dissimilar from a CPC scenario. The choice of parameters [2H1 is given in Table El 
(last column). In the definition of the CPX scenario ^S] the relations /i = 4msusY and 
|A,b| = I'^gl = 2?7isusY are fixed. Here, ■n^susY = 500 GeV is chosen. The parameter 
1712 is set to 200 GeV. Additionally the complex phases of At^h and rrig are fixed at 90° 
degrees. Variants of the CPX scenario are investigated to check the stability of the CPX 
results with respect to the choice of its parameters. The phases of A,b and rrig varied 
from from 0° to 180°, fi in between 500 and 4000 GeV. The scenario with arg(y4t,b) = 90° 
has very different features from a CPC case and therefore has good properties for a CPV 
benchmark scenario. 

The benchmark scan databases, containing masses, cross-sections and branching ratios 
for all three neutral Higgs bosons for a variety of different input parameters, are generated 
using both CPH 28j, a modified version of SUBHPOLE, and FEYNHIGGS 2.0 |46]. 
They are implemented in a modified version of HZHA ^^1- Initial-state radiation and 
interference between Higgsstrahlung and boson fusion processes are taken into account 
by HZHA. The parameter tan/3 is scanned from 0.6 to 40, and mH± is scanned from 4 to 
1000. In this region both Hi and H2 have a width below 1 GeV, negligible with respect 
to the experimental resolution of several GeV. 

Also in the CPV scenarios there is good agreement between the data and the back- 
ground estimation. The largest observed excess of the data over the background appears 
at rriHi = 40 GeV and mn^ = 105 GeV in the CPX benchmark set. The excess has 
(1 — CLfo) = 1.2 X 10~^, corresponding to a significance of 3.0 cr. Also here it should be 
noted that there is a large statistical probability of such an excess to appear somewhere 
in the parameter space under study. Limits on the MSSM parameter space are derived. 

Fig. 1221 shows the combined exclusion result for the CPX scenario with all phases 
equal to 90°, msusY = 500 GeV and /i = 2 TeV. Fig. 1221 (a) shows both the expected 
and observed 95% CL exclusion areas in the plane of mni and niu^- For heavy mna, 
Hi resembles the SM Higgs boson (almost completely CP-even) with very little effect 
from CP violation. The limit on the allowed mass of Hi for large mn^ is found to be 
mHi > 112 GeV. In the region below nin^ ^ 130 GeV CPV effects play a major role. 

Fig. I2S1 (b) shows the 95% CL exclusion areas in the parameter space of tan/5 and 
TTiHa- 0^6 can see that tan/5 < 2.8 is excluded. The band at tan/3 < 2.8 is excluded by 
searches for the SM-like Hi, while the band at tan/3 > 10 and mn^ < 120 GeV is excluded 
by searches for ZH2 and H1H2 topologies. 

Fig. 1221 (c) displays the parameter space of tan/? and rriHi. The range tan/? < 2.8 is 
excluded, and a lower limit of tan/3 > 3.2 exists if tt^Hi is below 112 GeV. For 4 < tan/3 < 
10, ZH2 production is dominant. The large difference between the expected and observed 
exclusion regions in the area of 4 < tan /3 < 10 is mainly due to a less than 2a excess in 
the data between m-h ~ 95 GeV and rrih ~ 110 GeV |7j, which corresponds to the mass of 
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H2 in this region. For ttihi < 50 GeV there are also unexcluded regions in the expected 
exclusion, which is due to dominant ZH2— >ZHiHi production with relatively large ttihi, 
yielding broad mass resolutions and therefore reduced sensitivity. 

In Fig. ESI (d) the exclusion area is shown in the parameter space of the theoretical 
input parameters tan/5 and r?7,H±, which are varied during the scan. Since the CPX 
scenario yields mn^ i=a ^h± for most of the scan points, this is very similar to Fig. EHl (b). 

The uncertainty inherent to the two theoretical approaches, CPH and FEYNHIGGS, 
is illustrated in parts (e) and (f ) of Fig. |2S1 The largest discrepancy occurs for large values 
of tan/?, where the FEYNHIGGS calculation (part (f)) predicts a higher cross-section for 
Higgsstrahlung, and hence a better search sensitivity than the CPH prediction (part(e)). 

The large impact of the value of the top quark mass on the exclusion limits is shown 
in Fig. 1^ For rritop = 179.3 GeV, the excluded range in tan/3 shrinks to tan/3 < 2.4. 

The effect of different choices of the CPV phases is illustrated in Figs. |2S1 and OBI 
Values of arg(y4t,b) = arg(mg) from 0° to 180° are displayed. Fig. |2H1 shows exclusion 
regions in the parameter space of tan/3 and mni for arg(74t,b) = arg(mg) = 90°, 60°, 30° 
and 0°. The lower limit on tan/3 > 1.9 in the scenario with phases of 60° is the lowest 
limit on tan/? in the CPV scenarios. At 30° and at 0° all areas for low ttihi and low 
tan /3 are excluded. The exclusion for the maximally CPV scenario CPX with 90° is very 
different from the exclusion of a CPC scenario (arg(y4t,b) = arg(??ig) = 0°). A variation of 
the second main parameter governing the size of CPV effects, "msusY, has similar effects 
on the exclusion to those of a variation of arg(y4t,b) = arg(mg). 

Fig. 1211 shows exclusion regions in the parameter space of tan /3 and mni for phases 
of (a) 135° and (b) 180°. The scenario in (a) is phenomenologically still similar to the 
original CPX scenario. The scenario in (b), which is in fact a CPC case, exhibits two 
allowed regions, of which the lower one from tan/3 = 3 to tan/3 = 13 has a low HiZ 
coupling. The unexcluded "hole" in the exclusion region for 90 < rriHi < 100 GeV is due 
to an excess of the background in the SM-like channels. 

Since the CPX scenario has a relatively high value oi fi = 2 TeV, which influences 
the mixing of the CP eigenstates into the mass eigenstates (see Eq. (jB))), fi is varied from 
Ij = 500 GeV to /i = 4 TeV in Fig. EH For /i = 500 GeV (Fig. EH (a)) and /i = 1 TeV 
(Fig. 123 (b)) the CPV effects are small. Therefore no unexcluded regions occur at small 
TT^Hi- The scenario with /i = 4 TeV (Fig. |2ZI (d)) has strong mixing and a suppression of 
pair production at large tan /?, resulting in an exclusion area that is considerably smaller 
than in the CPX scenario fFig. 1271 (c)). 

The proposal of the CPX scenario in J2Hj leaves the choice of msusY open, as long as 
the relations |v4t,b| = 2msusY, |"^g| = 2msusY and /i = 4msusY are preserved. In order to 
test the dependence on msusY, two scenarios are tested: Fig. |2H1 (a) shows the scenario 
CPXi.o, where the ratio between the parameters in the CPX proposal is preserved, while 
msusY is increased from 500 GeV to 1 TeV. Only small differences with respect to the 
CPX scenario with msusY = 500 GeV can be seen. Fig. |2H1 (b) shows the CPX scenario as 
given in Table HH but with only msusY set to 1 TeV, while the values of |A,b|, ling] and 
/i are kept flxed. This results in a decrease of the CPV effects and thus no unexcluded 
regions at small tt^Hi are observed. 
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8 Summary 

The searches for neutral Higgs bosons described in this paper are based on all data col- 
lected by the OPAL experiment, at energies in the vicinity of the Z resonance (LEPl 
phase) and between 130 and 209 GeV (LEP2 phase). The corresponding integrated lu- 
minosities are of about 720 pb~^. The searches addressing the Higgsstrahlung process 
e'^e~^TCZ and those for the pair production process e"'"e~— >-7ii?i2 are statistically com- 
bined. None of these searches reveals a significant excess of events beyond the predicted 
background level, which would indicate the production of Higgs bosons. 

From these results, model-independent limits are derived for the cross-section of a 
number of event topologies that could be associated to Higgs boson production. These 
limits cover a wide range of Higgs boson masses and are typically much lower than the 
largest cross-sections predicted by the MSSM. 

The search results are also used to test a number of "benchmark scenarios" of the 
MSSM, with and without the inclusion of CP-violating effects. 

In the CP-conserving case, new benchmark situations are investigated as compared 
to earlier publications. These are motivated either by new measurements of the b— i>s7 
branching ratio and the muon anomalous magnetic moment {g — 2)^, or in anticipation 
of the forthcoming searches at the proton-proton collider LHC In all these scenarios the 
searches conducted by OPAL exclude sizeable domains of the MSSM parameter space, 
even in those situations where the sensitivity of the LHC experiments is expected to be 
low. An overview of the results is given in Table E| In the "mh-max" scenario which, 
among all scenarios predicts the widest range of itly, values, the following limits can be 
set at the 95% confidence level: m\i > 84.5 GeV and ttia > 85.0 GeV; furthermore, if 
the top quark mass is fixed at the current experimental value of 174.3 GeV, the range 
0.7 < tan/3 < 1.9 GeV can be excluded (this range shrinks for higher values of rritop)- 

For the first time, a number of CP-violating MSSM scenarios are studied experimen- 
tally, where the CP-violating effects are introduced in the Higgs potential by radiative 
corrections. The "CPX" benchmark scenario is designed to maximize the phenomenolog- 
ical differences in the Higgs sector with respect to the CP-conserving scenarios. In this 
case the region tan/5 < 2.8 is excluded at 95% confidence level but no universal limit 
is obtained for either of the Higgs boson masses. However, for tan/5 < 3.3, the limit 
mHi > 112 GeV can be set for the mass of the lightest neutral Higgs boson of the model. 
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Luminosity table for e'''e -^7iZ production 


Channel Name 


Energies 


Luminosity 


Mass range 


Described 


nz^ 


(GeV) 


(pb-i) 


(GcV) 


in 


LEP 1 Channels 


qqrr/rrqq 


91.2 


46.3 


mn = - 70 


[nmni 


{nini^qqqq)ui/ 


91.2 


46.3 


m-H^ = 10 - 75, m-H^ =0-35 


mM 


qqvu 


91.2 


46.3 


mT^i = - 70 


[nmnj 


qqll 


91.2 


46.3 


mn = 2{)- 70 


12, ig 


LEP 2 Channels 


bbqq 


161-172 


20.4 


m^ = 40 - 80 


uniin] 


bbz^i^ 


161-172 


20.4 


m^^ = 50 - 70 


MM 


rrqq, bbrr 


161-172 


20.4 


m^^ = 30 - 95 


uniin] 


Xee 


161-172 


20.4 


m^ = 35 - 80 


mrni] 


X/i/z 


161-172 


20.4 


rriT^ = 35 - 80 


unmn 


bbqq 


183 


54.1 


m^Y = 40 - 95 


m 


(WiWi^4b)qq 


183 


54.1 


mH2 = 40 - 80, m-H^ = 10.5 - 38 


m 


hhi^i^/iniHi^ih)!^!^ 


183 


53.9 


m^^ = 50 - 95 


m 


bbrr/ 


183 


53.7 


m7^ = 30 - 100 


m 


rrqq/(?t:i?t:i^4b)rr 










bbee, bb^^u 


183 


55.9 


m7^ = 60 - 100 


m 


bbqq 


189 


172.1 


mn = ^0- 100 


'8 


(WiWi^4b)qq 


189 


172.1 


mn2 = 40 - 100, m-H^ = 10.5 - 48 


lai 


bbz^i//(WiWi^4q)z^i^ 


189 


171.4 


m-H = 50 - 100 


[H] 


bbrr/ 


189 


168.7 


m-H = 30 - 100 


m 


rrqq/(?t:i?t:i^4b)rr 










bbee, hhfifi 


189 


170.0 


m-H = 70 - 100 


m 


low mx{HiHi){vv ,ee,^^) 


189-192 


201.7 


mn^ = 45 - 90, mn^ = 2 - 10.5 


m 


bbqq 


192-209 


421.2 


m-^ = 80 - 120 





(HiHi^4b)qq 


192-209 


421.2 


mn2 = 80 - 120, mm = 12 - mn^/^ 


14.1.11 


bbi^i/ 
bbbbi^z^ 


192-209 
199-209 


419.9 
207.2 


m-H = 30- 120 
mn2 = 100 - 110, m-Hi = 12 - 771-^2/2 


7 


14.1.21 


bbrr / rrqq 


192-209 


417.4 


m^ = 80 - 120 


ill 


bbee, hh^ji 


192-209 


418.3 


m-^i = 40 - 120 





low mx{T~LiHi){vv ,ee,^i^jL) 


196-209 


396.9 


mn2 = 45 - 90, mn^ = 2 - 10.5 


m 




LEP 


2 Flavour-Independent Channels 




qqqq 


189 


174.1 


m^^ = 60 - 100 


isni 


qqz^i^ 


189 


171.8 


m-H = 30 - 100 


EOI 


qqrr, rrqq 


189 


168.7 


m-H = 30- 100 


isni 


qqee, qq^/x 


189 


170.0 


m-H = 70- 100 


3(1 


qqqq 


192-209 


424.2 


mw = 60 - 120 


(31l 


qqz^i^ 


192-209 


414.5 


m-K = 30 - 110 


ISH 


qqrr, rrqq 


192-209 


418.9 


m7^ = 60 - 115 


m 


qqee, qq^^u 


192-209 


422.0 


m-H = 60 - 120 


EH 



Table 1: List of the searches for the Higgsstrahlung process, 
reference or section where the search is described. 



The last column gives the 
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Luminosity table for 


e"'"e -^7ii7i2 production 




Channel Name 


Energies 


Luminosity 


Mass range 


Described 


7^17^2^ 


(GeV) 


(pb-i) 


(GeV) 


in 


LEP 1 Channels 


6b 


91.2 


27.6 


m-H^ = 40 — 70, m-}-i-^ =5 — 35 


.12,121 


qqrr, rrqq 


91.2 


46.3 


mn^ = 12 - 75, mn, =10-78 


[niinj 


6r, 4r2q, 2T4q 


91.2 


46.3 


m^a = 30 - 75, m-Hi =4-30 


naiini 


LEP 1.5 Channels 


4b 


130-136 


5.2 


S = 80-130,A = 0-50 


nn 


6b 


130-136 


5.2 


m7^2 = 55 — 65, rriT-ii > 27.5 


m 


LEP 2 Channels 


4b 


161 


10.0 


S = 80-130,A = 0-60 


uniiiii 


6b 


161 


10.0 


mn2 = 55 - 65, mm > 20.0 


unmu 


bbrr, rrbb 


161 


10.0 


m-H^ = 40 - 160, mm = 52 - 160 


MM 


4b 


172 


10.4 


S = 80-130,A = 0-60 


uniin] 


6b 


172 


10.4 


mm = 55 — 65, mm =25 — 35 


uniin] 


bbrr, rrbb 


172 


10.4 


mm = 37 - 160, mm = 28 - 160 


uniin] 


4b 


183 


54.1 


S = 80-150,A = 0-60 


M 


6b 


183 


54.1 


mm = 30 - 80, mm =12-40 


M 


bbrr, rrbb 


183 


53.7 


S = 70 - 170, A = - 70 


9' 


4b 


189 


172.1 


S = 80-180,A = 0-70 


m 


6b 


189 


172.1 


mm = 24 - 80, mm =12-40 


M 


bbrr, rrbb 


189 


168.7 


S = 70 - 190, A = - 90 


IHI 


4b 


192 


28.9 


S = 83-183,A = 0-70 


14.2.11 


4b 


196 


74.8 


S = 80-187,A = 0-70 


14.2.11 


4b 


200 


77.2 


S = 80-191,A = 0-70 


14.2.11 


4b 


202 


36.1 


S = 80-193,A = 0-70 


14.2.11 


high rrifi^ 4b 


199-209 


207.3 


S = 120 - 190, A = - 70 


14.2.11 


low my_^ 4b 


199-209 


207.3 


S = 100 - 140, A = 60 - 100 


14.2.21 


6b 


199-209 


207.3 


S = 90 - 200, A = 40 - 160 


14.2.31 


bbrr, rrbb 


192 


28.7 


S = 10- 174,A = 0-182 


14.2.41 


bbrr, rrbb 


196 


74.7 


S = 10- 182,A = 0- 191 


14.2.41 


bbrr, rrbb 


200 


74.8 


S = 10 -182,A = 0-191 


14.2.41 


bbrr, rrbb 


202 


35.4 


S = 10 -174,A = 0-182 


14.2.41 


bbrr, rrbb 


199-209 


203.6 


S = 70-190,A = 0-90 


14.2.41 



Table 2: List of the searches for pair production. The last column gives the reference or 



section where the search is described. The symbols S 
denote the Higgs mass sum and difference. 



mm + mm and A = m^m - mm 
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mn2 


mm 


Efficiency for the 


process 7i2Z^bbbbqc 


at yi 


(GeV) 


(GeV) 


192 GeV 


196 GeV 


200 GeV 


202 GeV 


206 GeV 


100. 


12. 


0.689 


0.684 


0.717 


0.733 


0.693 


100. 


20. 


0.651 


0.639 


0.653 


0.659 


0.586 


100. 


30. 


0.460 


0.461 


0.461 


0.470 


0.480 


100. 


40. 


0.270 


0.260 


0.283 


0.315 


0.323 


100. 


48. 


0.328 


0.325 


0.361 


0.392 


0.400 


105. 


12. 


0.538 


0.658 


0.702 


0.709 


0.701 


105. 


20. 


0.562 


0.618 


0.697 


0.658 


0.681 


105. 


30. 


0.490 


0.525 


0.509 


0.536 


0.497 


105. 


40. 


0.407 


0.306 


0.309 


0.316 


0.319 


105. 


50. 


0.433 


0.368 


0.355 


0.359 


0.370 


110. 


12. 






0.637 


0.682 


0.720 


110. 


20. 






0.625 


0.646 


0.532 


110. 


30. 






0.556 


0.549 


0.565 


110. 


40. 






0.380 


0.328 


0.343 


110. 


53. 






0.395 


0.341 


0.358 



Table 3: Efficiencies of the standard e"'"e^-^ZH-^qqbb analysis for the 

e"'"e~-^Z7Y2^Z7-^i7-^i— >qqbbbb final state (see Section r4.i.ij) . The uncertainties from 
Monte Carlo statistics are of the order of ±0.015. 



cut 


data 


tot. bkg. 


qq(7) 


tot 4-f 


Eff. 

7777^2=105 

mn,=20 


Eff. 
lH.2-^'HiHi 

777-^2 = 105 

mm =40 


Eff. 
'H2-^hh 

7777^2 = 105 


(l)-(5) 


503 


424.77 


123.94 


297.99 


0.75 


0.74 


0.59 


2-jet topology subsample A 


(6) 
(7) 

(8) 


371 
213 
135 


308.56 
201.77 
126.29 


108.66 

24.70 

22.85 


197.17 
177.07 
103.43 


0.75 
0.70 
0.68 


0.16 
0.15 
0.14 


0.53 
0.50 
0.49 


ANNa 


11 


10.0 


2.63 


7.39 


0.59 


0.11 


0.40 


4-jet topology subsample B 


(6) 


118 


112.32 


14.42 


97.80 


0.008 


0.57 


0.06 


ANNb 


8 


7.20 


2.83 


4.37 


0.003 


0.55 


0.05 


Total (2+4 jets) 


Sum 


19 


17.2 ± 0.6 


5.46 


11.8 


0.59 


0.66 


0.45 



Table 4: Cut Row in the missing energy analysis for 7i2i'iy—>-hhi'iJ and 
H2i^i^^'HiHii'i^^hhhhh'h' for 100 < 7777^2 ^ HO GeV for data taken at y/s = 199 to 
209 CeV (see Section U.l.'Jp . The uncertainty quoted for the total background is from 
Monte Carlo statistics only. 
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Efficiency 


A and B for 7i2^^^bbz/z/ decays 




ruH, (GeV) 


Efficiency of selection A 


Efficiency of selection 


B 


100. 


0.382 


0.044 




105. 


0.396 


0.050 




110. 


0.379 


0.050 




Efficiency A and B for l-L2Vv^flH.il-Livv^A:)\)bhvi' decays 


mn, (GeV) 


rriH, (GeV) 


Efficiency of selection A 


Efficiency of selection 


B 


100. 


12. 


0.686 


0.0 




100. 


20. 


0.561 


0.001 




100. 


30. 


0.280 


0.301 




100. 


40. 


0.090 


0.522 




100. 


48. 


0.195 


0.436 




105. 


12. 


0.707 


0.0 




105. 


20. 


0.587 


0.0 




105. 


30. 


0.349 


0.254 




105. 


40. 


0.113 


0.550 




105. 


50. 


0.179 


0.487 




110. 


12. 


0.677 


0.0 




110. 


20. 


0.585 


0.001 




110. 


30. 


0.402 


0.189 




110. 


40. 


0.109 


0.555 




110. 


50. 


0.131 


0.537 




110. 


53. 


0.186 


0.495 





Table 5: Efficiencies of the selections A and B of the missing energy analysis for 7-^2~^bbz/z/ 
and ?i2-^'^i'^i— ^bbbbz/z/ for 100 < m-}^^ < 110 GeV for data taken at y/s = 199 to 
209 GeV is used (see Section \'4.1.2Jj) . The uncertainty from Monte Carlo statistics is of the 
order of ±0.010. 
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Cut 


Data 


Total bkg. 


qq(7) 


4-fermi. 


Efficiency (%) 
m-H^ = m-H2 = 90 GeV 


l-L2l-Li^Ah Channel 


(1) 


39367 


39375.6 


30958.7 


8325.6 


99.8 


(2) 


13792 


13895.2 


8914.4 


4976.0 


98.4 


(3) 


4682 


4509.6 


1110.5 


3397.9 


88.3 


(4) 


3997 


3994.5 


707.7 


3285.7 


86.4 


(5) 


3474 


3431.0 


566.3 


2863.6 


85.6 


(6) 


3331 


3271.5 


520.4 


2749.9 


83.7 


£HiW2 > 0.95 


22 


19.9 ± 0.3 


6.5 


13.4 


49.4 


T^s^i^bbr+r- Channel 


Pre-sel 


336.0 


354.6 


96.0 


258.5 


53.8 


£rt2-tti > 0.64 


13 


13.2 ±0.4 


0.7 


12.4 


42.5 



Table 6: Cut Row in the 'H1H2 channels for high rriT-i-^ (see Section \4.2\) and for all data 
taken at y/s = 192 to 209 GeV: effect of the cuts on the data and the background, 
normalised to the integrated luminosity of the data. The two-photon background, not 
shown separately, is included in the total background. The signal efficiencies are given in 
the last column for m.-^^ =t^H2 =90 GeV. 





Efficiency for the process 7i2'^i^bbbb (high /tt-t^J at 
v/i = 206 GeV 


mn, (GeV) 


30.0 


40.0 


50.0 


60.0 


70.0 


80.0 


90.0 


mn, (GeV) 




30.0 


0.001 














40.0 


0.0008 


0.004 












50.0 




0.110 


0.215 










60.0 


0.103 




0.274 


0.364 








70.0 




0.254 




0.381 


0.388 






80.0 


0.200 




0.384 




0.425 


0.470 




90.0 




0.319 




0.388 




0.472 


0.479 


100.0 






0.374 




0.432 




0.435 


110.0 




0.341 




0.371 








120.0 


0.261 




0.349 




0.399 






130.0 




0.253 




0.368 








140.0 


0.231 




0.290 










150.0 




0.177 












160.0 


0.116 















Table 7: Efficiencies of the 7ii7i2— >bbbb analysis for high rra-^^ (see Section \4.2.1^ . The 
uncertainty from Monte Carlo statistics is typically of the order of ±0.015. The table is 
showing the Monte Carlo points produced. 
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Systematic uncertainties at ^/s- 


= 206 GeV 








W2Wi^bbbb (high rriT^J 


TC2'Hi^hhhh (low m-n-J 


W2Wi^bbr+r-,r+r-bb 


Source 


Signal eff. 


Background 


Signal eff. 


Background 


Signal eff. 


Background 


Detector modelling 


0.9% 


8.0% 


1.1% 


8.6% 


1.0% 


1.0% 


B-had. Decay Mult. 


0.9% 


1.2% 


2.1% 


1.9% 


1.0% 


3.0% 


B-had. Fragment. 


1.8% 


1.5% 


2.7% 


1.5% 


1.8% 


1.5% 


C-had. Fragment. 


0.0% 


0.5% 


0.0% 


0.0% 


0.0% 


0.0% 


4f-cross-section 


0.0% 


1.6% 


0.0% 


1.0% 


0.0% 


1.9% 


MC-Generators 


0.0% 


2.6% 


0.0% 


2.6% 


0.0% 


0.0% 


T identification 


0.0% 


0.0% 


0.0% 


0.0% 


0.9% 


15.0% 


LH modelling 


0.9% 


2.3% 


1.0% 


2.1% 


0.0% 


0.0% 


MC statistics 


2.0% 


5.0% 


3.0% 


5.2% 


1.0% 


3.0% 


Combined 


3.1% 


10.3% 


4.7% 


10.9% 


2.6% 


15.8% 



Table 8: Systematic uncertainties on the signal efEciency and background at -\/i= 
for the processes ?ii7i2— ^bbbb with high m-j-i^, low m-}^^ and for 7Y2'^^i— ^-bbr+r" 
(see Sections KTnKT^ and KT^ . 



206 GeV 
bb 



T^T' 



Cut 


Data 


Total bkg. 


qq(7) 


4-fernii. 


Efficiency (%) 
mn, = 30, mn^ = 100 GeV 


7ii7i2^4b Channel for low m^^^ 


(1) 


18519 


17802.0 


13705.2 


4096.9 


99.6 


(2) 


6538 


6427.8 


3971.6 


2456.3 


96.2 


(3) 


4215 


4048.0 


2082.4 


1965.6 


94.2 


(4) 


3618 


3497.5 


1546.8 


1950.6 


93.6 


(5) 


2712 


2625.9 


1188.5 


1437.4 


90.5 


(6) 


2477 


2389.4 


1060.0 


1329.5 


83.5 


£-«irt2 > 0.98 


8 


10.4 ± 0.1 


6.1 


4.3 


36.9 



Table 9: Cut Eow in the 7iiH2 channel for low mj^^ (see Section \4.2.2^) and for all data 
taken at y/s = 199 to 209 GeV: effect of the cuts on the data and the background, 
normalized to the integrated luminosity of the data. The two-photon background, not 
shown separately, is included in the total background. The signal efficiencies are given in 
the last column for m,-}i^=30 GeV and m,'}-i^=100 GeV in the 7^i7Y2^bbbb channel. 



Efficiency for the process 7i2'^i^bbbb (low m-^J at 
v/i = 206 GeV 


nin, (GeV) 


12.0 


20.0 


30.0 


m„, (GeV) 




90. 

95. 

100. 

105. 

110. 


0.269 
0.286 
0.305 
0.310 
0.298 


0.330 
0.341 
0.366 
0.358 
0.351 


0.370 
0.384 
0.369 
0.369 
0.366 



Table 10: Efficiencies of the 7ii7Y2— ^bbbb analysis for low m-}^-^ (see Section \4. 2. 2^) . The 
uncertainty from Monte Carlo statistics is of the order of ±0.011. 



33 





Efficiency for 7-^2'^i— ^bbbbbb at 
v/i = 206 GeV 


mn, (GeV) 


12.0 


20.0 


30.0 


40.0 


45.0 


50.0 


60.0 


70.0 


mn, (GeV) 




80. 


0.002 


0.188 


0.390 


0.465 










90. 


0.002 


0.263 


0.447 


0.595 


0.569 








100. 


0.001 


0.283 


0.486 


0.594 


0.639 


0.629 






no. 


0.001 


0.300 


0.552 


0.627 


0.662 


0.659 






120. 


0.002 


0.214 


0.512 


0.671 


0.664 


0.650 


0.695 




130. 


0.002 


0.292 


0.519 


0.635 


0.680 


0.670 


0.657 




140. 


0.000 


0.255 


0.536 


0.636 


0.646 


0.670 


0.649 


0.382 



Table 11: EfEciencies of the 7ii7i2— ^■'^^I'^^i'^^i-^bbbbbb analysis (see SectAon [4.2.!i^) . The 
uncertainty from Monte Carlo statistics is ±0.010. 





Efficiency for 7-^2'^i— ^bbr+r at 
v^ = 206 GeV 


mn, (GeV) 


30. 


40. 


50. 


60. 


70. 


80. 


90. 


mn, (GeV) 




30. 


0.0818 














40. 


0.1016 


0.2024 












50. 




0.3296 


0.3690 










60. 


0.3152 




0.4086 


0.3668 








70. 




0.3802 




0.3936 


0.4631 






80. 


0.4265 




0.3856 




0.4343 


0.4247 




90. 




0.3599 




0.4285 




0.3925 


0.4305 


100. 


0.3854 




0.4367 




0.3951 


0.4298 




110. 




0.3628 




0.4227 








120. 


0.3631 




0.3609 




0.4020 






130. 




0.2996 




0.3180 








140. 






0.3049 










150. 




0.2278 












160. 
















170. 

















Table 12: EfEciencies ofthe7ii7i2—>-hhT^T analysis (see Section \'4.2.4\) . The uncertainty 
from Monte Carlo statistics is of the order of ±0.02. 



Channel 


Total 
Background 


Data 


Section 


7^2Z-^bb(bb)qq 


135.4 


140 


14.1.11 


lH.2Vv-^hh{hh)vi' 


39.5 


36 


14.1.21 


7^2^i^bbbb(bb) 


19.9 


22 


14.2.11 14.2.21 I4.:^.3I 


Hirti^hhT+T- 


13.2 


13 


14.2.41 



Table 13: Typical numbers of background and data for the searches described in this 
publication. 
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Parameter 


no mixing 
/(no-mixing (2 TeV)) 


TTih-max 
/(mh-max"*") 
/(C mh-max) 


large-/! 


gluophobic 


small aoff 


CPX 


Parameters varied in the scan 


tan/3 


0.4-40 
/(0.7-40) 


0.4-40 


0.7-40 


0.4-40 


0.4-40 


0.6-40 


mA (GeV) 


0-1000 


0-1000 


0-400 


0-1000 


0-1000 


- 


mH± (GeV) 


- 


- 


- 


- 


- 


4-1000 


Fixed parameters 


mt (GeV) 


174.3 


174.3 


174.3 


174.3 


174.3 


174.3 


■msvSY (GeV) 


1000 
/(2000) 


1000 


400 


350 


800 


500 


M2 (GeV) 


200 


200 


400 


300 


500 


200 


f, (GeV) 


-200 

/(200) 


-200 

/(200) 
/(200) 


1000 


300 


2000 


2000 


mg (GeV) 


800 


800 


200 


500 


500 


1000 


Xt (GeV) 





VS^JT-SUSY 

/(Vemsusv) 
/(-\/6msusY) 


-300 


-750 


-1100 


At — fj. cot f3 


A,b (GeV) 


Xt + fi cot f3 


Xt -1- /i cot f3 


Xt + n cot /3 


Xt + fi cot /3 


Xt + ^i cot P 


1000 


arg(At,b) 

















90° 


arg(mg) 

















90° 



Table 14: Parameters of benchmark scenarios considered. Note that the values for Xt 
and y4t b are given for the MS-renormahzation scheme. For a description of the choice of 
parameters see Section Columns 2 to 6 refer to the CPC benchmark sets and the last 
column refers to the basic CPV benchmark set CPX. 



Limits on the MSSM scenarios 



Benchmark set 



Lower limit on mh (GeV) 



Lower limit on toa (GeV) 



Excluded tan/3 



no mixing 

no mixing (2 TeV) 

TTih— max 

TTih— max"*" 

constr. 77ih— max 

gluophobic 

small tteff 

CPX 



64.0 (60.0) 
83.3 (88.0) 
84.5 (88.5) 
84.5 (88.0) 
84.0 (88.0) 
82.0 (87.0) 
79.0 (83.0) 



84.3 (88.8) 
85.0 (89.0) 
84.0 (89.5) 
85.0 (89.0) 
87.5 (90.5) 
90.0 (95.0) 



0.8 < tan/?< 6.2 
0.9 < tan/3 < 4.2 
0.7< tan/3< 1.9 
0.7 < tan /3 < 1.9 
0.6 < tan/3 < 2.2 

tan/3<6.0 
0.4< tan/3< 3.6 

tan/3< 2.8 



(0.9 < tan/3 < 7.2) 
(0.9 < tan/3 < 4.3) 
(0.7 < tan/3 < 1.9) 
(0.7 < tan /3 < 1.9) 
(0.6 < tan/3 < 2.2) 
(tan/3 <8.0) 
(0.4 < tan/3 < 3.6) 
(tan/3 < 2.8) 



80.0 < mh < 107.0 
(81.0 < mh < 107.0) 



Allowed regions in the "large /i" scenario 



87.0 < niA 
(87.0 < toa) 



large ^ 



tan/3 > 6 
(tan/3 > 12) 



Table 15: Limits on m,\^ m,A and tan/? for the various benchmark sets. The median 
expected limits in an ensemble of SM background- only experiments are listed in paren- 
theses. The lower limits on rr^h and mp^ in the no mixing (2 TeV) scenario are only valid 
for rriA > 2 GeV. 
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Figure 1: Distributions of some ANN input variables for subsample A for 
e"'"e~-^7i2Z^bbz/z/ and e"'"e~-^7ii7iiZ^bbbbz/z/ (see Section \4.1.2^) in the region 
100 GeV < m-H^ < HO GeV for the data taken at ^/s = 199 - 209 GeV. Light grey 
(yellow) is the four-ferniion background, dark grey (green) the contribution from the 
qq(7) background. The dashed line shows the arbitrarily scaled signal expectation for 
m,fi^ = 12 GeV and m-^j = 105 GeV, as given in the legend. The points with error bars 
are the data. 
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Figure 2: Distributions of some ANN input variables for subsample B for 
e"'"e~— >-7i2Z^bbz/z/ and e"'"e~— ^T^iT^iZ— i>bbbbz/z/ (see Section \4.1.^) in the region 
100 GeV < m-H^ < 110 GeV. 
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Figure 3: Distributions of ANNp^ (a) and ANN-q (b) in the dedicated selection 
for e"'"e~— s>?i2Z— s>bbz/z/ and e'^ e~ —^HiT-CiZ^hhhhuh' (see Section \4.1.2\] in the region 
100 GeV < m-}^^ < 110 GeV. The arrows indicate the cut on the ANN output value. 
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Figure 4: Distributions of the reconstructed masses in the selection A and B of the missing 
energy channel for data taken at ^/s = 199 to 209 GeV in the region 100 GeV < m-)^,^ < 110 
GeV fsee Section WTl) . 
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Figure 5: Distributions of the likelihood input variables for the e'^e — i>7Yi7Y2— ^bbbb 
searches for high m-^-^ at 192-209 GeV (see Section [4.2.1]) . 
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Figure 6: Distributions of the input variables for the e'^e ^7Yi7Y2— ^bbbb searches for 
low m-H, at 199-209 GeV (see Section KTl) . 
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Figure 7: Searches for e+e~-^7ii7i2^bbbb with high m-^j (see Section \4.2.1[\4.2.'^) at 192- 
209 GeV. Likelihood outputs for a) high m-j^^^ , and b) low mj^^ . The arrow indicates the 
cut position. The signal is scaled with a factor of 50 with respect to a e"'"e~^7ii7i2— >bbbb 
signal for cos^(/3 — a) = 1. 
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Figure 8: Searches for e~^e —^l-iiH2 with high niy^-^ at 192-209 GeV. The sum of the 
reconstructed Higgs boson masses, m^' 



+ TTi^^, for (a) the 'H1H2 — >4b channel with the 



I rpr rpr I 



'Hi 



di-jet pairing combination which yields the smallest mass difference, AM = 
(b) the I-L1H2 —>-4b channel with the intermediate AM combination, (c) the I-L1H2 —>-4b 
channel with the largest AM combination. OPAL data are indicated by points with error 
bars, the four-fermion background by the light grey (yellow) histograms, and the two- 
fermion background by the darker grey (green) histograms. Shown as dashed histograms 
are the contributions expected from a Higgs boson signal with full strength at m^^j = 
m>^2 = 80 GeV for a luminosity of 207 pb~^ taken at ^/s = 199 — 209 GeV, added 
to the background expectation for a luminosity of 424 pb~^ and the cross-section for 
cos^(/3 — a) = 1. 
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Figure 9: Searches for e^e^^7iiH2 with low rriy^^ at 199-209 GeV. The sum of the 
reconstructed Higgs boson masses, nf-^^+Trf^^, for (a) the 7ii7i2 — ^4b channel with the di- 
jet pairing combination which yields the smallest mass difference, AM = |/;t^ 
the 7ii7i2 — *4b channel with the medium AM combination, (c) the 7ii7i2 - 
with the maximum AM combination. OPAL data are indicated by points with error bars, 
the four-fermion background by the light grey (yellow) histograms, and the two-fermion 
background by the darker grey (green) histograms. Shown as dashed histograms are the 
contributions expected from a Higgs boson signal with m?^^ = 30 GeV, m^-j-c^ = 100 GeV, 
added to the background expectation for a luminosity of 207 pb^^ and the cross-section 
for cos^(/9 — a) = 1. 
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Figure 10: Input variables j^ for the e"'"e -^7ii7i2^bbr+r searches at 192-209 GeV. 
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Figure 11: Likelihood output of the e'^e~^HiT-C2^hhT^r^ search at 192-209 GeV. The 
arrow indicates the cut position. The signal is scaled with a factor of 100 with respect to 
a e"'"e~— s>7^i7Y2^bbr+r^ signal for cos^(/? — a) = 1. 
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Figure 12: 
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normalization. The dashed line indicates the kinematic limit for y^ = 206 GeV. 
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Figure 13: Model-independent upper bounds on ax BR for (a) the e~^e —^7ii7i2^hhT^T 
channel and (b) the e"'"e~^7ii7i2^'^i'^i'^i— ^bbbbbb channel. The MSSM cross-section 
for T-iiH2 production with cos^(/3 — a) = 1 is taken as normalisation. The dashed line 
indicates the kinematic limit for ^/s = 206 GeV. 
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Figure 14: Results for the no mixing benchmark scenario. The figure shows the excluded 
regions in darker grey (green) and theoretically inaccessible regions in light grey (yellow) 
as functions of the MSSM parameters in four projections: (a) the (m,^, m^) plane, (b) the 
fmA, tan/3j plane, (c) the (m.^, tan/3j plane and (d) the (mii±, tan/5j plane. The dashed 
lines indicate the boundaries of the regions expected to be excluded at the 95% CL if only 
SM background processes are present. The region excluded by Yukawa searches, Z-width 
constraints or decay independent searches is shown in dark grey (red). In (b) the hatched 
area is still allowed by BR(b— >S7) searches. 
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Figure 15: Results for the no mixing (2 TeV) benchmark scenario described in the text 
of Section \7.1l See Fig. \T^ for the notation. 
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Figure 16: Results for the mi^—raax benchmark scenario. See Fig. ^| for the notation. 
The dotted hne in (c) shows the observed 99.9% confidence limit. The differently shaded 
regions in (b) show the exclusion for different values of mtop, as written in the Plot. The 
upper limit on m\^ for different values of mtop is shown in (c), as expressed in the plot. 
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Figure 17: Results for the mh— max"*" benchmark scenario described in the text of Sec- 
tion \7.1l See Fig. E| for the notation. This scenario is excluded by the BR(b-^s'~f) 
constraint for ttia < 600 GeV. 
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Figure 18: Results for the constrained mY^—raax benchmark scenario described in the text 
of Section \7.1\ See Fig. E| for the notation. The differently shaded regions in (b) show 
the exclusion for different values ofrritop, as written in the Plot. The upper limit on rrih 
for different values ofmtop is shown in (c), as also expressed in the plot. 
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Figure 19: Results for the large /x benchmark scenario described in the text of Section 1 7. 11 
See Fig. ^^for the notation. 
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Figure 20: Results for the gluophobic benchmark scenario described in the text of Sec- 
tion \7.1\ See Fis.[T^for the notation. The hatched area in (c) is allowed by the BRlb^sj) 
constraint. 
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Figure 21: Results for the small Oefr benchmark scenario described in the text of Sec- 
tion \7.1\ See Fig. \T^for the notation. 
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Figure 22: Diagram illustrating the effective coupling of a Higgs mass eigenstate to the 
Z in Higgsstrahlung. The complete mass eigenstate Tii is composed of admixtures of 
h, H and A. Here the h, H and A denote the CP-even and CP-odd weak eigenstates, 
respectively. Only the CP-even admixtures h and H couple to the Z, while the CP-odd A 
does not. Therefore the coupling of the mass eigenstate is reduced with respect to a CPC 
scenario. 
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Figure 23: The CPX MSSM 95 % CL exclusion areas. Excluded regions are shown for 
(a) the ('mHi,mH2J plane, (b) the (mii^,ta.ia(3) plane, (c) the frnHj ,tan /3j plane and (d) 
the (m}i± ,tan f3) plane. Figure (e) shows the fmHi ,tan f3) of the CPH calculation alone, 
(f) shows the same projection of the FEYNHIGGS 2.0 calculation. See Fig. El for the 
notation. The dash-dotted line in (c) shows the area excluded on the 99.9% confidence 
level. In (b) and (d) the area excluded by Z width constraints or by decay independent 
searches is too small to be displayed. 
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Figure 24: The CPX MSSM 95 % CL exclusion areas in the (mi{i,ta.nj3) plane, using scans 
with (a) nit = 179.3 GeV and (b) rrit = 169.3 GeV. Due to the change in the top masses 
a strong difference is observed compared to Fig. \2^ (c). See Fie:. flU for the notation. 
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Figure 25: The CPX MSSM 95 % CL exclusion areas in the (mii^,ta.nP) plane, using 
scans with (a) arg(/lt,b) = ai'g('^g) = 90°, (b) arg(y4t,b) = arg(mg) = 60°, (c) arg(74t,b) = 
arg(mg) = 30° and (d) arg(y4t,b) = arg(^g) = 0°. While the CPV phases decrease, effects 
from CP violation like the strong H2^>HiHi contribution vanish. See Fig. El for the 
notation. 
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Figure 26: The CPX MSSM 95 % CL exclusion areas in the (m}ii,ta.nP) plane, using scans 
with (a) a.Tg{At^h) = arg(mg) = 135° and (b) arg(y4t,b) = arg(mg) = 180°. See Fie:. [Til for 
the notation. 
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Figure 27: The CPX MSSM 95 % CL exclusion areas in the (mn^jtanP) plane, using 
scans with (a) /i = 500 GeV, (b) /i = 1000 GeV, (c) /i = 2000 GeV (CPX) and (d) 
/i = 4000 GeV. See Fig.\Mfor the notation. 
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Figure 28: The CPX MSSM 95 % CL exclusion areas in the fmHutan/Jj plane, using 
scans (a) preserving the CPX ratios of fi, A^^t and msusY, using mgusY = 1 TeV, and (b) 
'^susY = 1 TeV keeping fj, and Ab^t at their CPX values. See Fig. \T^ for the notation. 



62 



